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A candidate of a neutron-halo nucleus, 31 Ne, contains a single neutron in the pf 
shell. Within the Glauber and eikonal models, we analyze reactions used to study 
31 Ne. We show in a 30 Ne+n model that the magnitudes of the total reaction and 
above all of the one-neutron removal cross sections of 31 Ne on 12 C and 208 Pb targets 
strongly depend on the orbital angular momentum of the neutron, thereby providing 
us with efficient ways to determine both the spin-parity and structure of the ground 
state of 31 Ne. Besides these inclusive observables, we also calculate energy and 
parallel-momentum distributions for the breakup of 31 Ne, and show their strong 
dependence upon the orbital of the valence neutron in the bound state of 31 Ne. 
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I. INTRODUCTION 

Exploring nuclei near the neutron and proton driplines is making rapid progresses in 
and beyond the p, s<i-shell region. The Ne isotopes raise interesting structure problems. 
The alpha cluster structure around 20 Ne is known for many years [l|, 0]. Recently 17 Ne, 
an 15 0+p+p Borromean system, has been found to have a large charge radius due to a 
significant amount of s 2 component 0]. For the very neutron-rich Ne, Na, and Mg isotopes 
with N ~ 20, one of the most important issues is the vanishing of the shell gap, which 
causes a mixing of normal and intruder configurations, and has significant influence on the 
properties of those nuclei ji-lQ]. The importance of deformation around 30 Ne is stressed in 
Refs. 0, @|, in contrast to the result of a mean- field calculation 11]. The heaviest Ne isotope 
synthesized so far is 34 Ne. It may be a dripline nucleus considering that 33 Ne is unstable to 
neutron decay [l2T ]. 
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The nucleus 31 Ne with N = 21 neutrons attracts our special attention in view of its pos- 
sible halo structure containing a lps/2 and/or O/7/2 valence neutron. Its neutron separation 
energy S n is 0.33 MeV, though it has large uncertainty [l3j]. The ground state spin-parity 
of 31 Ne is thus expected to be either 3/2" or 7/2". The former possibility may happen 
because the single-particle energy of the neutron orbit with low orbital angular momentum 
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receives a considerable shift near the neutron dripline [14|, [15 
a shell model |I| and one using a microscopic cluster model of 
inversion. 

The rare isotope 31 Ne was first produced in a projectile fragmentation reaction |11 
Nowadays, an intense beam provided by the Radioactive Ion Beam Factory (RIBF) at 
RIKEN can produce 31 Ne in sufficiently large amounts (several particles per second). Very 
recently, the total reaction cross sections or of heavy Ne isotopes on 12 C target [l8j and the 
one-neutron removal cross sections cr_ n of 31 Ne on 12 C and 208 Pb targets have been measured 
for the first time around 230 MeV/nucleon The purpose of this article is to analyze 
the sensitivity of <tr, <r_ n and other dissociation cross sections to the orbit of the 31 Ne 
valence neutron. During the completion of this theoretical work, the one-neutron removal 
cross sections of 31 Ne measured at RIKEN became available 19[. We seize this opportunity 



to compare our calculations with the data to draw conclusions about the structure of the 
ground state of 31 Ne. 

We describe 31 Ne as a system consisting of a 30 Ne core (c) and a weakly-bound valence 
neutron (n). The core is assumed to be in its + ground state though its excitation energy is 
fairly low. Considering that structure model we evaluate the total reaction and one-neutron 



20H23| on both light ( 12 C) and heavy 



removal cross sections within the Glauber formalism 
( 208 Pb) targets, and compare the values obtained for the 1^3/2 and O/7/2 possible configu- 
rations of the 31 Ne ground state. To predict the sensitivity of more exclusive observables 
(e.g. energy and parallel-momentum distributions) to the ground state configuration, we 



also perform calculations within the eikonal model |20H22l. |24j . Since both light and heavy 



targets are considered, we use the Coulomb correction to the eikonal model (CCE) j25-27l 

This article is structured as follows: After a summary of the Glauber and eikonal for- 
malisms (Sec. [II]), we detail the inputs of our calculations in Sec. IIHI Our results and analysis 
are presented in Sec. OV] Sec. [VTJ contains the conclusions and perspectives of this study. 



II. THEORETICAL FRAMEWORK 



As mentioned in the introduction, we consider in this study two reaction models. First, 
the Glauber model 20|-|23| is used to evaluate the total reaction and one-neutron removal 
cross sections of 31 Ne. Second, the eikonal model 
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221 . |24j is used to compute the disso- 
ciation cross section as a function of the energy and parallel-momentum between the 30 Ne 
core and the neutron after breakup [ 
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Both models are based on Glauber's idea |24| 
to describe the influence of the collision onto the initial projectile-target wave function by a 
multiplying amplitude e lx , 



(1) 



where the phase x is assumed to depend only on the transverse component b of the projectile- 
target relative coordinate. In the present work, this phase is obtained by folding a profile 
function that describes nucleon-nucleon effective interactions with the projectile and target 
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densities. In the eikonal approximation, however, it is more usual to derive it from optical 
potentials that simulate the interaction between the projectile constituents and the target. 
In a general interpretation of the eikonal model [20|, the adiabatic approximation em- 



ployed in the Glauber model is not assumed, which invalidates the simple ansatz ([T]) [30 
The adiabatic approximation ignores the excitation energy compared to the incident energy, 
leading to a well-known unphysical result for the Coulomb dissociation. In order to solve 
this problem and still maintain Eq. ([I]), we only need to correct the Coulomb phase appro- 



priately |25|, |26| . This approximate version is the CCE and its accuracy has been tested by 
comparison to the exact eikonal calculation in Ref. [27| . 

In this section, we briefly present both approaches, emphasizing their common points and 
differences that make them complementary. 



A. Glauber formalism 



Provided that 31 Ne can be seen as a neutron loosely bound to a 30 Ne core whose wave 
function is the same as that of an isolated 30 Ne, <7_ n can be obtained from the difference 



between the projectile and the core interaction cross sections |20l . 131 



a_ n ( 31 Ne) = ai( 31 Ne) - a^Ne). 
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Computing the interaction cross sections is not easy because it excludes inelastic scattering, 
which cannot be properly treated if no description of the internal structure of the projectile 
is considered. Fortunately, if the number of bound excited states is small, cxi can be well 
approximated by the reaction cross section ctr, which can be easily computed within the 
Glauber formalism j20|, 21, 23]. For 31 Ne, which has only one known bound state, i.e. its 
ground state, this approximation is legitimate. For 30 Ne, however, o"R( 30 Ne) will overestimate 
<Ti( 30 Ne) by cr inc i( 30 Ne), in which the projectile is excited towards its 2 + and 4 + bound 
excited states. Nevertheless, the inelastic scattering being a phenomenon occurring near 
the nuclear surface, its contribution is not expected to be significant at incident energies of 
200-300 MeV/nucleon where the surface transparency becomes large. The approximation 



<x_ n ( 31 Ne) « a R ( 31 Ne) - cr R ( 30 Ne) 



(3) 



seems thus reasonable. 

The Glauber model expresses the nuclear part of the reaction cross section for a nucleus 
X impinging on a target T as the integral of the reaction probability with respect to the 
transverse components b of the X-T relative coordinate 20, 21, 2^ 



OR 



)db, 



(4) 



where the phase-shift function \ models the nuclear interactions between the colliding nuclei. 
As mentioned earlier, in the Glauber formalism, this phase is expressed as a function of the 
densities of the target p T and the impinging nucleus px- It also depends on profile functions 
Tnn describing effective nuclear interactions between the nucleons. At the optical limit 
approximation of the Glauber model (OLA) the nuclear phase-shift functions are usually 



given by [20j, [21|, |23j 



x N (b) 



PT(r') Px (r")T NN (b - *' + s")dr"dr', 



(5) 
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where s' and s" are the transverse components of the internal coordinate of the target 
(r') and the impinging nucleus (r"), respectively. The OLA is therefore equivalent to the 
double-folding of an effective nucleon-nucleon interaction. Note that the profile functions 
Tnn depend on the nucleons considered: Their expression for identical nucleons (pp or nn) 
is not the same as for the proton-neutron (pn) interaction. Therefore, in our calculations, 
expression (jSJ) is actually split into 4 terms. This is done as follows: Replace Px^nn with 
Px-^Np + Px-^Nn using the proton and neutron densities of the projectile X and change px 
by pJp + Pt followed by renaming N of Tnn in accordance with the density. 

As shown by Abu-Ibrahim and Suzuki, the OLA (jSj) misses some higher-order terms, 
which can be included using the symmetrized expression 33, [34 



X N (b) = ^ ( 



dr"p x (r")T NN (b + s'-s" 
+ J dr"p x (r") jl - exp - J dr' p T (r')T NN (b - s' + a") X\ . 



(6) 



The Glauber calculations presented in the following are performed using this expression (jSJ) 
in Eq. (|J]). Again the actual phase-shift function in our calculations is split into four terms. 
The details about the calculation of the densities and the profile functions are summarized 
in Sec. HI 

For the carbon target, the Coulomb contribution to the total reaction cross section is 
neglected, the reaction being fully nuclear dominated. However, this may no longer be done 
for heavy targets. For the lead target, we add incoherently to the nuclear reaction cross 
section (T4|) the Coulomb contribution at first-order (see Sec. HVBl) . 



B. Coulomb-corrected eikonal description of reactions 

Since we are also interested in the influence of the 31 Ne structure on other observables, 
like energy and parallel-momentum distributions, we perform calculations within the eikonal 



model [20|, |2l|, |24j . Indeed this model enables us to compute differential cross sections con 



sidering both Coulomb and nuclear interactions, as their interferences [28|, |29|]. The eikonal 
model assumes a cluster structure of the projectile and usually describes the interaction 
between the clusters and the target by optical potentials. 

In this work, the projectile P ( 31 Ne) is assumed to be made up of a neutron n of mass 
m n initially bound to a 30 Ne core c of mass m c and charge Z c e. This two-body projectile 
is impinging on a target T of mass and charge Z^e. The neutron has spin / = 1/2, 
while both core and target are assumed to be of spin zero. These three bodies are seen as 
structureless particles. Fig. [TJ schematizes the set of coordinates we use in the following. 
The c-n relative coordinate is denoted by r and P-T relative coordinate by R, with Z and 
b its longitudinal and transverse components, respectively. In Fig. [TJ the transverse parts of 
the c-T (b c ) and n-T (b n ) coordinates are shown as well. 

The structure of the projectile is described by the internal Hamiltonian 

H = ^ + V cn (r), (7) 

where p is the relative momentum of the neutron to the core, p cn = m c m n /mp is the reduced 
mass of the core-neutron pair (with mp = m c + m n ), and V cn is the potential describing the 
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FIG. 1: Jacobi set of coordinates used for the eikonal calculations: r is the projectile internal 
coordinate, and R = b + ZZ is the target-projectile coordinate. The transverse components of the 
core-target and neutron-target coordinates are denoted by b c and b n , respectively. 



core- neutron interaction. This potential includes a central part, and a spin-orbit coupling 
term (see Sec. Illlfl . 

In partial wave Ij, the eigenstates of H are defined by 

H 0( f> ljm (E,r) = E<f> ljm (E,r), (8) 

where E is the energy of the c-n relative motion, and j is the total angular momentum 
resulting from the coupling of the orbital momentum I with the neutron spin I. The negative- 
energy solutions of Eq. ([8]) correspond either to the physical bound state of the projectile, or 
to orbitals occupied by the neutrons of the core, which are forbidden to the valence neutron 
by the Pauli principle. The former is denoted by (f>i o j omo (E ) in the following. These wave 
functions are normed to unity. The positive-energy states describe the broken-up projectile. 
Their radial part r^ 1 Ukij are normalized according to 

Ukij(r) — > cos SijFi(kr) + sm SijGi(kr), (9) 

r— >oo 

where k = ^2[i cn E /ft 2 is the wave number, 5ij is the phase shift at energy E, and Fi and 
Gi are respectively the regular and irregular Coulomb functions [Hj]. 

At the eikonal approximation, the amplitude appearing in Eq. (CQ) can be divided into 
three factors (27| 

e tx = e lx ^e lx e lx , (10) 

where the dependence on the transverse coordinate b has been omitted for clarity. The 
elastic Coulomb phase Xpt describes the projectile-target Rutherford scattering. It reads 



24 



X c PT (b) = 2 V ln(Kb), (11) 

where K is the wave number of the projectile-target relative motion and r] = 
ZtZ c 6 2 / (4neohv) is the P-T Sommerfeld parameter, with v the initial P-T relative velocity. 
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Besides the deflection of the projectile trajectory, the Coulomb interaction also con- 
tributes to the breakup of the projectile. Acting only on the core, it indeed induces a tidal 
force between both components of the projectile. The Coulomb phase X C in Eq. (JTOJ sim- 
ulates that tidal force (see e.g. Eqs. (16) and (17) of Ref. The slow decrease of this 
phase at large b leads to divergence in the calculation of the breakup cross sections j27j. To 
overcome this problem, Margueron, Bonaccorso, and Brink have proposed a correction to 
this Coulomb term [IBJ . It consists in replacing at first order the Coulomb phase x° by the 
first order of the perturbation theory x F ° ( see Eq. (22) of Ref. (27|) following 



->• e 



>x 



iX C + iX F °- 



(12) 



Because at large b the first-order phase x FO decays exponentially, correction flTJJ) solves the 
aforementioned divergence problem. In addition, it restores most of the missing dynamical 
effects in the eikonal model, which enables us to describe reactions taking on (nearly) the 
same footing both Coulomb and nuclear interactions at all orders (27| . 

In the eikonal model, the nuclear interactions between the projectile constituents and 
the target are usually described by optical potentials chosen in the literature. In that case, 
the nuclear phase x N is expressed as integrals over Z of these potentials (20I . 21, 24]. In 



the present case, no experimental data exist to constrain such a potential for the interaction 



between the Ne core and the target. Following Ref. 23], we approximate the nuclear phase 



for each projectile constituent by the OLA (j5J). Therefore 

x N (b,s) = x? T (b c ) + xZAK), 



(13) 



where XcT anc ^ XnT are respectively the c-T and n-T nuclear phases. They are computed 
using Eq. (j5J), in which the density px is replaced by the 30 Ne density or a Dirac delta 
function, respectively. 

To evaluate elastic-breakup cross sections within the CCE we proceed as explained in 
Ref. 27j. The elastic-breakup amplitude reads 



Hjm 



e ix -iX° + iX F °)e lx ' 



Hojomo (Eq 



(14) 



where 07 is the Coulomb phase shift [3| 

In the following, we consider two breakup observables. The first is the breakup cross 
section as a function of the c-n relative energy E after dissociation 
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The second breakup observable is the parallel-momentum distribution 



(15) 
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where 9 k = &rccos(k/ k\\) is the colatitude of the c-n relative wavevector k after breakup. 



(16) 
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100 


240 


1000 




0.65 


43.71 


72.53 


3.36 


2.58 


3.12 


1.54 


1.29 


1.38 


1.20 


0.70 


43.81 


73.55 


3.40 


2.59 


3.16 


1.56 


1.31 


1.39 




0.75 


43.85 


74.52 


3.46 


2.59 


3.20 


1.59 


1.33 


1.41 




0.65 


41.01 


68.83 


3.44 


2.74 


3.22 


1.58 


1.33 


1.41 


1.25 


0.70 


41.15 


69.79 


3.48 


2.73 


3.25 


1.60 


1.34 


1.43 




0.75 


41.22 


70.72 


3.52 


2.72 


3.28 


1.62 


1.36 


1.44 



TABLE I: Single-particle potentials for 30 Ne, rms radii of 30 Ne and total reaction cross sections of 
30 Ne+ 12 C collision at the incident energy of 100, 240 and 1000 MeV/nucleon. Lengths, energies, 
and cross sections are given in units of fm, MeV, and b. 

III. DENSITIES AND POTENTIALS 

The calculation of the cross sections described in the previous section requires projectile 



and target densities and profile functions. In our study, we follow Ref. [23 . 

We first construct 30 Ne densities. We assume the internal wave function of this nucleus 
to be a Slater determinant of single-particle orbitals generated from the following potential 

U(r) = -V f(r) + V,rl I • a~f(r) + V c (r)±^, (17) 

where T3 has eigenvalue 1 for neutrons and —1 for protons, and / is the Woods-Saxon form 
factor 

f(r) = {l + exp[(r-R)/a]}- 1 , (18) 

where radius R = tqA 1 ^ 3 , with A = 30. The spin-orbit strength is set to follow the system- 
atics 36l . 



V 1 = 22-U[(N-Z)/A]t 3 (19) 

in MeV. The Coulomb potential Vc is taken from a uniform charge distribution. The values 
of ro and a are varied around standard values, and Vq is determined separately for neutrons 
and protons to fit S n and S p . The resulting values are denoted V n and V P , respectively. 
The neutron and proton densities of 30 Ne, p° and pj?, are calculated from the occupied 



orbits by removing approximately the effect of the center of mass motion [23j. The root mean 
square (rms) radii for neutron, proton, and matter distributions (r", r£, r™) are listed in 
Tablefl] The table also contains cr R ( 30 Ne) for a 12 C target at 100, 240 and 1000 MeV/nucleon. 
The second value of the incident energy is chosen because it is close to that of the RIKEN 



experiment [18|, [19] , and that profile functions are available at that energy [37J . The choice 
of T NN is explained later in this section. 

Since 31 Ne is assumed to exhibit a 30 Ne-n cluster structure, its densities are obtained 
from the 30 Ne densities computed above, and the wave function 4>ijm f° r the 30 Ne-n relative 
motion. The latter is determined by solving the Schrodinger equation (jSJ) in either the lpz/2 
or O/7/2 orbital. The 30 Ne-n interaction is simulated by the same mean- field potential as for 
30 Ne ffTT|) . but with a different central depth Vq. 
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FIG. 2: (Color online) Single-particle energies of the (full lines) and O/7/2 (dashed lines) 

valence neutrons of 31 Ne as a function of Vb for three different diffuseness parameters a (fm). The 
ro value is 1.25 fm. Horizontal dotted-line at —0.33 MeV denotes the experimental energy. 



Fig. [2] displays the single-particle energies of IP3/2 and O/7/2, e{p) and e(f), as a function 
of Vo for three choices of diffuseness parameter a, the radius parameter being fixed to tq = 
1.25 fm. With increasing a, e(p) decreases very rapidly, whereas e(f) shows a mild change. It 
is therefore possible to obtain the expected shell inversion by considering a sufficiently large 
diffuseness (e.g. a = 0.75 fm). For actual calculations, the strength Vq is set to reproduce 
the S n value of 0.33 MeV (see Table HT|) . Note that these potentials are also used as V cn in 
the calculations of the wave functions 0/ jm that appear in the eikonal model (see Sec. Ill Bp . 

The proton and neutron densities of 31 Ne, p v p and pp, are calculated including the recoil 
effect, which means that the difference between the centers of mass of 31 Ne and 30 Ne is 
treated properly 

p P (r') = pP(r') (20) 
p P (r') = p»(r') + p n (r'), (21) 

where r' is the internal coordinate of 31 Ne. In these expressions, and p" are the contribu- 
tions of the 30 Ne core to the 31 Ne densities. They slightly differ from the densities of 30 Ne, 
because of the recoil effect 

Pl ,n {r') = / P^i^r + r'MUr^dr, (22) 

m 

where r is the 30 Ne-n relative coordinate. In Eq. fT2T]) p n denotes the contribution of the 
valence neutron to the 31 Ne density 

Pn(r') = E / 6 ^i r ~ r'Mi im (r)\ 2 dr. (23) 

m 

Fig. E] displays the matter density of 31 Ne (pp = p v p + p P ) as well as its contributions from 
the 30 Ne core (p c = ~p v c + p°) and the valence neutron (p n ). The IP3/2 orbit (left panel) 
reaches far in distances and extends the tail of the 30 Ne density significantly beyond 6 fm. 
On the contrary, the O/7/2 orbit (right panel) changes the density only slightly even near the 
surface. 
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FIG. 3: (Color online) Matter density of 31 Ne (full lines) and contributions from the 30 Ne core 
(dashed lines) and the valence neutron (dotted lines). The ro value is 1.25 fm, and a is 0.75 fm for 
the lps/2 orbit (left panel) and 0.70 fm for the O/7/2 orbit (right panel). 



Table HT1 lists the valence-neutron single-particle energies (e), the rms radii of 31 Ne for 
the neutron, proton and matter distributions (rf,, 7p, r™), and <jR( 31 Ne) for a 12 C target at 
100, 240 and 1000 MeV/nucleon. We also give the rms radius of the valence-neutron orbit 
r n = a/ (r 2 )- This r n value turns out to be around 7fm for the 1^3/2 orbit but, due to the 
larger centrifugal barrier, is much smaller for the Of 7/2 orbit: about only 4 fm. Interestingly, 
although the matter radii of 31 Ne and 30 Ne depend on the potential sets (see Tables [J and 
HT|) . their difference remains unchanged: Ar = r™ — r™ is 0.19-0.20 for 1^3/2 and 0.04 fm for 
O/7/2. The constancy of Ar within the set of the same / suggests that Ar is insensitive to the 
shape of the potential but determined by S n and I. Despite the fact that the single-particle 
energy is only — 0.33MeV, Ar is not very large even for because the mass number of 
the core nucleus is fairly large. 

Since the neutron separation energy S n of 31 Ne is not accurately known, we also perform 
calculations with a slightly deeper potential (see last line of Table [TT]) in order to examine 
the S n dependence of crR( 31 Ne) and cr_ n ( 31 Ne) values. This potential gives e(p) = —0.6 MeV 
instead of —0.33 MeV. The matter radius is reduced by only 0.05 fm, but the r n value 
changes by about 0.8 fm. The decrease of <jR( 31 Ne), and thus of <r_ n ( 31 Ne), on carbon is 
only about 10 mb. However, the cr_ n ( 31 Ne) value on lead is expected to be considerably 
reduced. We will discuss this in Sec. IIVBI 

The target densities used in our calculations are obtained from experimental data. For 
both 12 C and 208 Pb, the proton densities are derived from empirical charge densities by 
removing the finite size effect of protons. The neutron density of 12 C is obtained as explained 



in Ref. [23|. For 208 Pb, the neutron density is obtained by subtracting the proton density 
from the matter density [32| taken from a Hartree-Fock calculation. 

Other key inputs to compute the cross sections of Sec. [TT1 are the profile functions Tnn 
that correspond to effective nucleon-nucleon interactions. These functions are parametrized 



in the usual way |23|, [37 

WM = ^i<;exp(-^H, (24) 
47rp NN \ zp NN / 

where ct^n is the total cross section for the N-N collision, onn is the ratio of the real to 
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ro a V e(p) e(/) r n rf> r p p rf ctr 

100 240 1000 



0.65 52.28 -0.22 -0.33 4.25 3.40 2.59 3.16 1.58 1.32 1.41 

1.20 0.70 51.82 -0.33 -0.05 7.20 3.66 2.60 3.35 1.70 1.41 1.50 

0.75 50.87 -0.33 - 7.35 3.72 2.60 3.40 1.72 1.43 1.52 

0.65 47.65 -0.01 -0.33 4.39 3.48 2.74 3.26 1.61 1.35 1.43 

1.25 0.70 47.87 -0.24 -0.33 4.47 3.52 2.73 3.29 1.63 1.37 1.45 

0.75 47.41 -0.33 -0.07 7.44 3.79 2.73 3.48 1.75 1.45 1.55 

1.25 0.75 48.52 -0.60 -0.55 6.62 3.71 2.73 3.43 1.74 1.44 1.54 



TABLE II: Properties of the potentials describing 31 Ne. Using various potential geometries, we 
adjust S n = 0.33 MeV in either the 1^3/2 orbit or the 0f 7 / 2 one. Last-line potential reproduces 
S n = 0.60 MeV in the 1^3/2 orbital. Rms radii of the corresponding densities are listed as well as the 
total reaction cross sections of 31 Ne+ 12 C at incident energies of 100, 240, and 1000 MeV/nucleon. 
Lengths, energies, and cross sections are given in units of fm, MeV, and b. 

the imaginary part of the N-N scattering amplitude, and /?nn is the slope parameter of 
the N-N elastic differential cross section. The values of these parameters are taken from 



Ref. [37J. Note that they differ for the interaction between identical nucleons (pp or nn) and 
for the interaction between a proton and a neutron (pn). To analyze the sensitivity of our 
calculations to this choice of profile functions, we also perform calculations that ignore the 
difference between pp (or nn) and pn interactions. In those tests, we use the parameters of 
Tnn given in Ref. [231 ] . 

The profile functions (124]) combined to the densities of 30,31 Ne and of the target enable us 
to compute the phase-shifts ([6]) for the Glauber calculation. The same parameters are used 
to derive the OLA fl5]) used in the Coulomb-corrected eikonal calculation. To this end, the 
densities of the projectile and the target are expanded on a Gaussian basis 

p( r ) ~ Ci exp ( ~ \ air2 ) ■ ( 25 ) 

i 

This enables us to solve analytically the integrals appearing in Eq. ([5]) and partly in Eq. ([6]). 
The values q and are available from the authors. 

In the eikonal model, the nuclear phase ((5]) is added to the elastic Coulomb phase ( TTTT) 
and the corrected Coulomb phase (jl~2]) to obtain the eikonal phase (jl~0|) . That phase is then 
numerically expanded into multipoles of rank A. To this end, we use a Gauss quadrature 
on the unit sphere similar to the one considered to solve the time-dependent Schrodinger 
equation in Ref. {38|. The number of points along the colatitude is set to Ng = 12, and the 
number of points along the azimuthal angle is A^ = 30 in most cases but goes up to 40 
when large As are considered. For the carbon target, the calculation requires a rather large 
number of multipoles: A max = 16 in the l]?3/2 case, and A max = 12 in the O/7/2 one. For 
the lead target, a smaller number of multipoles is needed: A max = 8 for the 1^3/2 state, and 
Amax = 6 for the 0f 7/2 one. 

The eigenfunctions of the projectile Hamiltonian Hq ([8]) are computed numerically with 
the Numerov method using 1000 radial points equally spaced from r = up to r = 100 fm. 
This rather large value is required in order to reach convergence in the radial integrals 
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FIG. 4: (Color online) Total reaction cross section of 31 Ne on a 12 C target as a function of incident 
energy. Dotted lines are the results with the OLA phase shifts ([5]) . See caption of Fig. [3] for ro and 
a. 

appearing in Eq. (fl~4"l) and in the calculation of r n , the rms radius of the valence neutron (see 
Table IITI) . The integrals over b appearing in Eqs. (fT5l) and (fl~6l) are performed numerically 
from b = up to b = 400 fm with a step Ab = 1 fm. In the 1^3/2 case this integral had to 
be done up to 600 fm to reach convergence when a lead target was considered. 




IV. DISCUSSION OF THE lp 3/2 AND 0/ 7/2 ASSUMPTIONS WITHIN THE 

GLAUBER MODEL 

A. Total reaction cross sections 

Fig. H] compares <7R( 31 Ne) on a 12 C target calculated within the Glauber model (see 
Sec. Ill Aj) for the IP3/2 (full line) and O/7/2 (dashed line) orbits as a function of the 31 Ne 
incident energy. The phase-shift function is calculated using Eq. fl6]). The projectile density 
is obtained using the potential sets of radius ro = 1.25 fm with diffuseness a = 0.75 fm for 
the 1^3/2 orbit and a = 0.70 fm for the O/7/2 orbit. At all energies the relative difference in 
<7r between both configurations is about 5-10%. For example, as listed in Table HTl <7R( 31 Ne) 
at 240 MeV/nucleon is 1.45 b for the p orbit and 1.37 b for the / one. Thus the difference 
of <jR( 31 Ne) depending on whether the orbital angular momentum of the valence neutron 
is 1 or 3 amounts to 87 mb. Though not very large, this difference may be sufficient to 
determine which assignment is favorable in comparison with experiment [181 ] . 

The reaction cross section is larger for a 1^3/2 neutron than for the O/7/2 neutron because 
the integral appearing in the phase shifts ([6]) extends on a larger domain in the former case 
than in the latter. This variation in <jR( 31 Ne) with the projectile configuration, being mostly 
due to the change in the valence-neutron orbital is therefore rather small. Indeed, most of 
<jR( 31 Ne) is contributed by the 30 Ne core, whose reaction cross section does not vary much 
with the potential set: <jR( 30 Ne) = 1.36 b for r = 1.25 fm and a = 0.75 fm, and 1.34 b for 
ro = 1.25 fm and a = 0.70 fm. On the contrary, the increase in <tr due to the addition of the 
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FIG. 5: (Color online) Nuclear contribution to the total reaction cross section of 31 Ne on a 208 Pb 
target as a function of incident e nerg y. The dotted lines are the results obtained with the average 
profile functions taken from Ref. [23|]. See caption of Fig. [3] for ro and a. 



valence neutron, erR( 31 Ne) — cTR,( 30 Ne), is strongly dependent on the assumed configuration: 
The increase turns out to be 96 mb for the lpz/2 orbit and 26 mb for the O/7/2 orbit at 
240 MeV/nucleon. Following Eq. ([3]), this result suggests the one-neutron removal cross 
section to be an observable more sensitive to the projectile configuration (see Sec. HVBj) . 

To investigate the sensitivity of our calculations to the construction of the phase-shift 
function, we also compute <TR,( 31 Ne) using the OLA (JSj) (dotted lines in Fig. H]). As is usually 



observed |23l. I37I]. the OLA tends to predict larger cross sections. However, the difference 
between the reaction cross sections obtained with the lpz/2 configuration and the O/7/2 one 
is about the same using OLA (J5|) as when the phase-shift function flU]) is used. 

Fig. displays <7R( 31 Ne) on a 208 Pb target calculated with only the nuclear phase shifts. 
The effect of Coulomb breakup is discussed in the next subsection. As observed for the car- 
bon target, the difference between the lpz/2 (full line) and O/7/2 (dashed line) configurations 
is small though non- negligible. As mentioned earlier, this difference comes mainly from the 
valence-neutron contribution. The increase of <tr from 30 Ne to 31 Ne is even more striking 
for a 208 Pb target. It is almost ten times larger considering a lpz/2 valence neutron than a 
O/7/2 one. At 240 MeV/nucleon, the reaction cross section increases from 4.36 b to 4.69 b 
in the former case while it goes from 4.33 b to only 4.37 b in the latter. 

Since the proton and neutron densities of the lead target are different, we examine how 
much the cross sections depend on the choice of the profile function T NN . Fig. [5] compares 
two sets of calculations, one which employs different interactions between pp (or nn) and 
pn (full and dashed lines), and the other which uses the averaged interaction taken from 
Ref. [23J (dotted lines). As observed in Fig. [5], the choice of the averaged interaction tends 
to slightly overestimate the cross sections below 300 MeV/nucleon. 

The enhanced cross section for the lp^/2 orbit reflects the spatial extension of the neutron 
orbit. If its S n value is increased to, say 0.6 MeV as shown in Table HTl o"R( 31 Ne) gets smaller 
compared to that with S n = 0.33 MeV: At 240 MeV/nucleon, it is reduced by 12 mb for 
carbon and by 65 mb for lead. These cross sections are still significantly larger than those 
for the O/7/2 neutron case. 
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FIG. 6: (Color online) One-neutron removal cross section of 31 Ne, approximated by <7R,( 31 Ne) — 
CR,( 30 Ne), on a 12 C target as a function of incident energy: (range between full lines) and 

O/7/2 (range between dashed lines). The experimental point is from Ref. [19|. 



B. One-neutron removal cross sections 

As mentioned in Sec. Ill we evaluate the one-neutron removal cross section cr_ n for 31 Ne 
on carbon and lead targets using approximation ([3]). Fig. [6] shows the results obtained on 
a 12 C target as a function of the 31 Ne incident energy for both (full lines) and O/7/2 

(dashed lines) configurations. To evaluate the sensitivity of these results to the potential 
set used to generate the projectile densities, we have performed the calculations with the 
different potentials given in Tables HI and HT1 Though the 1P3/2 or O/7/2 orbits vary with the 
potential set, they predict very similar cr_ n values: In both cases these values are contained 
between the pairs of lines shown in Fig. |6j Hereafter we use the potential set with ro = 1.25 
fm and a = 0.75 fm for the IP3/2 orbit, and the set with ro = 1.25 fm and a = 0.70 fm for 
the O/7/2 orbit unless otherwise mentioned. 

As discussed in the previous subsection, the interesting result of this set of calculations 
is that <r_ n is always much larger for a 1^3/2 valence neutron than for a O/7/2 one. At 



240 MeV/nucleon, close to the energy of the RIKEN experiment [19|, the former configu- 
ration leads to a cross section of about 96 mb, whereas the latter gives only 26 mb. This 
difference is basically due to the larger spatial extension of the p orbit compared to that of 
the / orbit, which is due to the change in the centrifugal barrier. The experimental cross 
section amounts to 79(7) mb [19:]. This value, being both close to our 1^3/2 calculation 
and much higher than our O/7/2 one, favors a ground state wave function for 31 Ne strongly 
dominated by a configuration in which the valence neutron is in the orbital coupled to 
a 30 Ne core in its + ground state. This comparison therefore suggests a 3/2~ spin-parity 
for the 31 Ne ground state, rather than the 7/2~ deduced from the naive shell model. 

As shown in Fig. EJ the difference in the magnitude of cr_ n increases at lower incident 
energies. An experiment performed at such an energy (e.g. a few tens of MeV/nucleon) 
would improve the confidence in the identification of the 31 Ne configuration. 

To evaluate <r_ n for a 208 Pb target we may no longer neglect the Coulomb contribution to 
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FIG. 7: Electric dipole transition strength of Ne as a function of the c-n relative energy E after 
dissociation. For the 1^3/2 transition, both contributions from the s- and c?-wave continuum states 
are shown. Results obtained with distorted waves (DW, full lines), plane waves (PW, dashed lines), 
and orthogonalized plane waves (OPW, dotted lines) are shown separately. 



the one-neutron removal process. Since the Coulomb interaction contributes mostly to the 
elastic breakup, we add an estimate of the Coulomb-breakup cross section to the reaction 
cross section computed within the Glauber framework. To this end, we use the first-order of 
the perturbation theory, considering only the dominant dipole transition. In that approxi- 
mation, the l]?3/2 neutron is excited to continuum states with / = or 2, whereas the O/7/2 
neutron is moved to d or g positive-energy states. This Coulomb contribution to <r_ n can be 
estimated by integrating the electric dipole transition strength dB(El)/dE multiplied by the 
photon number spectrum over the excitation energy j20[. Fig. [7] compares the dB(El)/dE 
distributions for the initial p (thick full lines) and / orbits (thin full line). In the former 
case, the partial-wave contributions are shown as well. Note that the result obtained from 
the initial O/7/2 configuration is multiplied by 10 for readability. These quantities depend 
on the choice of a minimum impact parameter 6 min from which the Coulomb breakup is 
assumed to contribute. However, the dependence of cr_ n on b min is found to be moderate 
around b mm = 12.7 fm, which is obtained from b m - m = r e ff (31 1 / 3 + 208 1 / 3 ), with r c g = 1.4 fm. 

The dipole strength obtained for the lps/2 configuration is concentrated at low excitation 
energy. The s wave gives a larger contribution to that distribution than the d wave at 
E < 0.5 MeV, but the d wave dominates over the s wave with increasing energy. On the 
contrary, dB(El)/dE for the O/7/2 initial state, besides being much smaller than the 
one, has a completely different energy dependence: It is flat and extends to high energies. 
This suggests that differential observables, like energy or parallel-momentum distributions, 
could be used to discriminate between these two possible configurations (see Sec. IV)) . 

To evaluate the sensitivity of this calculation to the c-n final state interactions, we evaluate 
the dipole strength for the initial IP3/2 bound state using distorted waves (DW, i.e., positive- 
energy eigenstates of the c-n Hamiltonian (jSJ); full lines), plane waves (PW; dashed lines), or 
orthogonalized plane waves (OPW, i.e., plane waves orthogonalized to the Pauli- forbidden 
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bound states of Hamiltonian (jSJ) [26[; dotted lines). Interestingly only the s wave contribution 
is sensitive to the continuum description: That value is much reduced in the vicinity of its 
maximum when DW are considered instead of PW or OPW. Nevertheless, these changes do 
not affect the results as much as to modify our conclusions. 

At 240 MeV/nucleon, and using DW, we obtain 0.81 b for the Coulomb contribution to 
cr_ n : 0.32 b from the s wave and 0.49 b from the d waves. This value is added incoherently 
to the nuclear contribution to cr_ n , which is estimated to be about 0.33 b in the Glauber 
model. The resulting cx_ n value turns out to be 1.14 b. As expected, the dipole strength 
obtained for the / orbit is much smaller: Its contribution to o"_ n is a mere 57 mb. The 
nuclear contribution is evaluated in the Glauber model to be about 34 mb, leading to a total 
(j_ n = 91 mb. This is about one order of magnitude smaller than the cross section for the 



p orbit. The experiment performed at RIKEN gave o"_ n = 712(65) mb [19J. Thus again 
slightly below our theoretical prediction for the lps/2 configuration, and much higher than 
the cross section obtained for the OiV/a orbit. This confirms the shell inversion predicted 
by former structure calculations [H, 16], in agreement with the analysis of Nakamura et al. 
19|. Note that evaluations of the Coulomb contribution using PW or OPW lead to similar 



results: large <r_ n for the lp^/2 configuration, and small <r_ n for the O/7/2 one. 

As mentioned in the last paragraph of the previous subsection, the Coulomb breakup 
contribution will be very sensitive to S Q of the 1^3/2 orbit. We have repeated the calculation 
assuming S n = 0.6 MeV. The a_ n value for S a = 0.6 MeV is predicted to be 0.75 b, of which 
0.49 b is due to the Coulomb breakup. Changing S n from 0.33 MeV to 0.6 MeV thus reduces 
o"_ n by 0.32 b. This is much larger than the corresponding reduction (65 mb) in the nuclear 
breakup contribution. Since a_ n changes significantly as a function of S n mainly because of 
the Coulomb dissociation, a close analysis of a_ n on a 208 Pb target can give some constraint 
on the S n value of 31 Ne. The one-neutron removal cross section obtained with S n = 0.6 MeV 
being closer to the experimental value, suggests that the one-neutron separation energy of 
31 Ne might be higher than 0.33 MeV. However, this reduction from theory to experiment 
may also be due to a spectroscopic factor for the 1^3/2 configuration smaller than one. Other 
observables, like energy or parallel-momentum distributions for elastic breakup, may provide 
further valuable information. 



V. EIKONAL CALCULATION OF ENERGY AND PARALLEL-MOMENTUM 

DISTRIBUTIONS 

Besides the significant difference in magnitude between the one-neutron removal cross 
section, the distinction between the 1^3/2 and O/7/2 configurations for 31 Ne could be made by 
looking at differential breakup observables, like energy or parallel-momentum distributions. 
To analyze the influence of the 31 Ne configuration on such cross sections, we perform elastic- 
breakup calculations within the Coulomb-corrected eikonal model (CCE, see Sec. Ill Bl and 
Ref. [271]). Unlike the Glauber model, the CCE solves the divergence problem posed by the 
Coulomb interaction between the projectile and the target. This enables us to take account 
of nuclear and Coulomb interactions on the same footing and to include their interference in 
the description of the reaction process. The following calculations are performed with the 
inputs detailed in Sec. IIHI 

The elastic-breakup cross sections obtained for 31 Ne impinging on a carbon target at 
240 MeV/nucleon are shown in Fig.[S]as a function of the energy E between the 30 Ne core and 
the neutron after dissociation. The total cross section for the IP3/2 configuration is displayed 
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FIG. 8: Energy distribution for the elastic breakup of 31 Ne on carbon at 240 MeV/nucleon. The 
total cross section for the initial IP3/2 bound state is plotted as well as its major partial- wave 
contributions. The result obtained with the initial O/7/2 bound state, multiplied by 10, is shown 
for comparison. 



with the thick full line, while its dominant s—g contributions are plotted with interrupted 
lines. The breakup cross section obtained considering the O/7/2 ground state is depicted 
with the thin full line. Note that it is multiplied by 10 for readability. Both distributions 
differ significantly. First, as already mentioned in Sec. IIVBI the magnitude of the O/7/2 cross 
section is much lower than the 1^3/2 one. Second, the 1^3/2 distribution is strongly peaked 
at low energy, whereas the O/7/2 distribution extends over a broader energy domain. This 
confirms that in addition to one-neutron removal cross sections, energy distributions could 
be used to determine the configuration of 31 Ne ground state. 

The two bumps observed in the O/7/2 cross section at about 5 and 9 MeV correspond 
to /5/2 and g 9 / 2 resonances of widths r / 5/2 — 1.5 MeV and r 0gg/2 — 3 MeV, respectively. 
These resonances are produced by the c-n potential used in this calculation (see Table ITTl) , 
but were not fitted to any known state. In the present work they have thus no physical 
meaning. However, this result indicates that if 31 Ne were to exhibit resonant states with a 
strong 30 Ne-n cluster structure, these could be revealed by a measurement of elastic breakup 



on a light target [39 



These resonances are also present in the 1^3/2 calculation, but the bumps they generate 
are less marked than in the O/7/2 case. The 1^3/2 orbit, being two quanta of orbital angular 
momentum further away from the resonances than the O/7/2 state, is indeed less prone to 
be excited towards that part of the continuum. 

We also perform a similar calculation for a 208 Pb target. The corresponding energy 
distributions are plotted in Fig. |9j As in the nuclear breakup case, the two configurations lead 
to very different results. Not only is the magnitude of the distribution strongly dependent 
on the initial state (note that the O/7/2 cross section is multiplied by 10), but also its shape 
clearly reveals the configuration of 31 Ne ground state. As in Figs. |E]and[7J the 1^3/2 energy 
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FIG. 9: Same as Fig. [8] for a lead target. Note the change in the energy axis. 
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FIG. 10: Parallel-momentum distribution for the elastic breakup of 31 Ne on carbon at 
240 MeV/nucleon. The cross section for the initial IP3/2 bound state is compared to that ob- 
tained with the initial O/7/2 bound state. The latter is multiplied by 10. 



distribution is peaked at low energy and decreases rapidly with E. The 0/ 7 / 2 cross section, 
on the contrary, is much natter. 

Another observable that is often used to discriminate the orbital of valence nucleons is the 
parallel- momentum distribution 40N42| . In that case, the breakup cross section is evaluated 



as a function of the parallel-momentum between the core and the neutron after dissociation. 
Fig. [10] depicts the parallel- momentum distribution for the elastic breakup of 31 Ne on a 
carbon target at 240 MeV/nucleon. The results obtained with both the 1^3/2 (thick line) 
and O/7/2 (thin line) configurations are shown. Note that here also the latter is multiplied 
by 10 for clarity. 

The signature of the initial configuration is even clearer here than in the energy distri- 
bution. Besides the significant change in magnitude, we observe that the 0/ 7 / 2 parallel- 
momentum distribution is much broader than that of the 1^3/2 configuration. This distribu- 
tion can be understood as a reminiscence of the initial bound-state wave function expressed 



in the momentum space |40|, |42|. The large spatial expansion of the 1^3/2 wave function 



translates into a narrow momentum distribution, which is revealed in this breakup cross 
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section. On the contrary, the narrower spatial distribution of the O/7/2 state leads to the 
broader parallel-momentum distribution observed in Fig. [101 



VI. CONCLUSION AND PERSPECTIVES 

The very neutron-rich isotope 31 Ne (N = 21) is located in a region where mixing of normal 
and intruder shell configurations is expected. In a naive shell model, the 31 Ne ground state 
would thus be seen as a 30 Ne core in its + ground state to which a O/7/2 valence neutron 
is added. However, some calculations predict this valence neutron to be in a 1^3/2 intruder 
orbital instead [H, lf|. If this were the case, the low angular momentum of the orbital 



combined to the low one-neutron separation energy of 31 Ne (S n ~ 0.33 MeV [13]) would 
suggest this nucleus to exhibit a one-neutron halo. 

Recently, the new RIBF facility at RIKEN has produced a 31 Ne beam at about 
230 MeV/nucleon. This beam is sufficiently intense to allow the measurement of its to- 



tal reaction and one-neutron removal cross sections on carbon and lead targets [18|, [19] . The 
present work aims at analyzing the sensitivity of these cross sections to the structure of the 
exotic isotope 31 Ne. To this aim we use the Glauber model detailed in Ref. |23] to evaluate 
ctr and cr_ n . This theoretical work shows that both ctr and cr_ n computed considering a 
IP3/2 configuration for 31 Ne are larger than those obtained with a O/7/2 valence neutron. Es- 
pecially, the difference in <r_ n is significant enough to doubtlessly discriminate between the 
two possible configurations. During the completion of this theoretical work, the one-neutron 
removal cross sections of 31 Ne measured at RIKEN became available The comparison 
of these data to our calculations suggests a strong lpz/2 configuration in the wave function 
of 31 Ne ground state, confirming, independently from the analysis of Nakamura et al. 19], 
the expected shell inversion in 31 Ne. We therefore conclude the spin-parity of that ground 
state to be 3/2~ rather than 7/2" as suggested by the naive shell model. 

Since other observables could be used to test this shell inversion, we have also performed 



prospective calculations within the Coulomb-corrected eikonal approximation [27| for the 
breakup of 31 Ne on both carbon and lead targets. These calculations confirm that a O/7/2 
configuration would lead to much smaller breakup cross sections than if the valence neu- 
tron were in the intruder lps/2 orbital. They also show that the shape of the energy and 
parallel-momentum distributions could be used to distinguish between the two possible con- 
figurations. Indeed, whereas assuming a IP3/2 valence neutron gives energy distributions 
peaked at low energy, the O/7/2 configuration leads to distributions that reach much higher 
energies. We have also observed that the parallel-momentum distribution is much narrower 
when the bound state is assumed in the p partial wave than in the / one. The measurement 
of these distributions would therefore provide a complimentary way to confirm the structure 
information obtained from the recent RIKEN measurement of er_ n . 

At such a distance from the valley of stability and near the region of the island of inver- 
sion, the 31 Ne ground state may not be composed of a single configuration. An extension 
of the reaction models used in this work to a multiple-configuration description of the pro- 
jectile structure, as the one proposed by Summers et al. (43], would definitely improve the 
reaction model. Such a model would indeed help understanding the influence of a multiple- 
configuration structure of the projectile upon reaction observables. 
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